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Abstract 

A gauged (2+l)-dimensional version of the Skyrme model is investigated. The gauge group is U(l) 
and the dynamics of the associated gauge potential is governed by a Maxwell term. In this model 
there are topologically stable soliton solutions carrying magnetic flux which is not topologically 
quantized. The properties of rotationally symmetric solitons of degree one and two are discussed in 
detail. It is shown that the electric field for such solutions is necessarily zero. The solitons' shape, 
mass and magnetic flux depend on the U(l) coupling constant, and this dependence is studied 
numerically from very weak to very strong coupling. 

PACS number(s): ll.10.Kk, ll.10.Lm, 11.27.+d, 12.39.Dc 



I Introduction 

The Skyrme model is a generalized non-linear sigma model in (3+1) dimensions ||T|. It 
has soliton solutions which, after suitable quantization, are models for physical nucleons 
0]. The theory is invariant under the group SO(3) lso of iso-rotations, and electromagnetism 
is introduced into the model by gauging a U(l) subgroup of SO(3) iso , see for details. 
The resulting fully coupled Skyrme-Maxwell system is mathematically hard to analyze, but 
of considerable physical interest: it is here that one should compute the Skyrme model's 
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prediction for the proton-neutron mass difference for example. In fact such a computation 
(which necessarily involves quantum theory) was first attempted in where the authors 
made various approximations based on the smallness of the fine structure constant. The 
result was a mass difference of m p — m n = 1.08 MeV, which has the wrong sign (in nature 
the proton is 1.29 MeV lighter than the neutron). 

In this paper we investigate classical properties of a (2+l)-dimensional version of the 
gauged Skyrme model. The model, to be introduced in section [H| is a gauged version of the 
baby Skyrme model studied in || and || and contains a dynamical abelian gauge field. It 
has soliton solutions which are stable for topological reasons and which carry magnetic flux. 
However, the gauge symmetry is unbroken and the solitons differ from the much studied flux 
tubes or vortices in the abelian Higgs model in that their magnetic flux is not quantized. 
Contrary to the situation in (3+l)-dimensional Skyrme-Maxwell theory it is possible to 
compute certain soliton solutions in our model explicitly with moderate numerical effort, and 
to investigate their structure quantitatively. Thus we study the dependence of the magnetic 
flux and the solitons' mass on the electromagnetic coupling constant, assess the back- reaction 
of the electromagnetic field on the matter fields and make some semi-quantitative statements 
about the long range inter-soliton forces. 

II Gauging the Baby Skyrme Model 

The model we want to study is defined on (2+1) -dimensional Minkowski space, whose sig- 
nature we take to be (—,+,+). Points in Minkowski space are written as (t, x) or simply 
x, with coordinates x a , a = 0, 1,2, and the velocity of light is set to 1. We will mostly be 
concerned with static fields and sometimes use the label i = 1,2 for the coordinates of the 
spatial vector x. The basic fields are a scalar field describing matter and a U(l) gauge 
potential A a for the electromagnetic field. More precisely, 4>(x) is a 3-component vector 
satisfying the constraint (f) ■ 4> = <p\ + 0| + 0i = 1 , thus lying on a 2-sphere which we denote 
by Si 

The model is a gauged version of the baby Skyrme model studied in detail in || and ||. 
The Lagrangian considered there is invariant under global iso-rotations of the field cf) about 
a fixed axis n e «S|. Taking n = (0, 0, 1) for definiteness such a rotation can be written in 
terms of the rotation angle x [0, 27r) as 

(0i, 02, h) -> (cos x 0i + sin % 02, - sin xh + cos x 02, 03) • (1) 

We write SO(2) iso for the group of such rotations. 

Here we couple electromagnetism to the baby Skyrme model by gauging the SO(2) iso 
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symmetry. Thus we require invariance under local rotations 

<t>^0{x)<t>, (2) 

where 0(x) is an S0(2) iso rotation matrix which depends on x. For infinitesimal rotation 
angles e(x), this becomes 

<fi — ► <fi + en x 0. (3) 

The abelian gauge field A a transforms as A a — > A a — d a e, so defining the covariant derivative 
via 

D a cf) = d a (f) + A Q n x 0, (4) 

one has D a (0(x)<p) = 0(x)D a (f>, as desired. Finally we define the curvature or field strength 
F a p = d a Ap—df3A a with electric components Ei = F i0 and the magnetic component B = F 12 . 
Thus we can write down the Lagrangian of our model: 

L = -H J d 2 x (^(D a (f>) 2 + ^ (D Q cP x Df,*f>) 2 + ^ 2 (1 - n • 0) + JLf%A . (5) 

The first term is a gauged version of the 0(3) sigma model Lagrangian (see [0]), the second 
is a gauged Skyrme term, the third term may physically be thought of as a pion mass term 
and the last term is the usual Maxwell Lagrangian. There are 4 free parameters in this 
model. H has the dimension energy, A and 1/fj, are of dimension length while g represents 
the coupling strength to the gauge field and is also of dimension length. We will discuss our 
choice of parameters in further detail below, but for the time being we fix our energy scale 
by setting H = 1. 

It is worth recalling that the Skyrme term is necessary in the (ungauged) Skyrme model 
to prevent soliton solutions from collapsing to singular spikes. In the gauged model, however, 
the Maxwell term has the same scaling behaviour as the Skyrme term, which suggests that 
there could be stable solitons in a "Skyrme-Maxwell theory without a Skyrme term" . This 
possibility is studied in 0. Here we retain the Skyrme term because we are also interested 
in properties of soliton solutions in the limit of vanishing electromagnetic coupling. The 
Skyrme term ensures the existence of stable solitons in this limit. 

The Euler-Lagrange equations for this model can be written conveniently in terms of 

J a = <p x D a (j> + X 2 D p <f)(D p <f) ■ 4> x D a $) (6) 
and the conserved current 

j a = n ■ J a . (7) 

They read: 

D a J a = |i 2 nx<^ (8) 
d a F«P = g 2 f. (9) 
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The second equation has three components which we want to write explicitly in terms of the 
electric and magnetic field, for later use. The (3 — component is Gauss's law: 



diEi=j . (10) 

The remaining two equations are particularly simple when expressed in terms of polar coordi- 
nates (r, 6) for x. Defining polar and radial coordinates of the current j via jg = X\ji — x 2 ji, 
jr = (xiji + x 2 j 2 )/r and analogously for the electric field E i: we obtain 



dE r 1 dB 2 



dt r d6 



= 93 



r 



dE e dB 2 . 
The energy E of a configuration (0, A a ) is the sum of the kinetic energy 

T = \J<P* ((^o^) 2 + y W x A0) 2 + j 2 E^ (12) 
and the potential energy 

V = ^Jd 2 x {{D^f + {D 2 <t>f + A 2 (A0 x D 2 <f>) 2 + 2/i 2 (l - n • </>) + Ifi 2 ^) . (13) 

In this paper we are only interested in finite energy configurations, so we require that for 
all t 

lim d>(t, x) = n . (14) 

r— »oo 

This boundary condition allows the Euclidean space R 2 to be compactified to a topological 
2-sphere iS 2 so that, at a given time t, fields <fi may topologically be regarded as maps 

<t> ■ Sl^Sl (15) 

Such maps are topologically classified by their degree Q, which is an integer and can be 
calculated via 

Q=^-f d 2 x<t>-(d 1 <t>xd 2 <t)) . (16) 

47T JS^ 

The degree is a homotopy invariant of <fi and therefore conserved during time evolution. 



Ill Static Solutions 

It is a well-known and important feature of the 0(3) sigma model and its generalization to 
baby Skyrme models that the potential energy of a configuration is bounded below by the 
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(2) 



(3) 



modulus of its degree (or a suitable multiple thereof). A similar result holds in our model, 
but its proof requires a little work. In ||, where the the potential energy functional 

M = \j d 2 x ((D^) 2 + (D 2 <(>) 2 + (1 - n ■ 4>f + B 2 ) (1) 

is studied in detail, it is shown that K ux [0, A] is bounded below by 47r|Q|. Thus, changing 
variables in the expression for V via x — > /xx, and discarding the positive definite Skyrme 
term we have the inequality 

V[<j>, M>\j d 2 x ^0) 2 + (D 2 4>) 2 + 2(1 - n -<t>) + ^B 2 ^ . 

Now note that, since < (1 — n •</>)< 2, it follows that (1 — n • </>) > |(1 — n • <fi) 2 . Thus 
we also deduce 

V[<f>, Ai]>\f d 2 x [{D^f + (D 2 <(>) 2 + (1 - n • 4>) 2 + ^B 2 ^ . 
If n/g > 1 it then follows immediately that 

V[4>,Ai] >K U3E [0,A] >4tt|Q|. (4) 
If /Vg — 1 011 the other hand we have 

2 2 

> ^Ku,[0,A] >4tt^|Q|. (5) 

In both cases we have therefore found a topological lower bound for the potential energy V. 

Our next goal is to find static configurations of given degree n > which minimize the 
potential energy V. To find these we exploit the symmetries of our model. Both V and Q are 
invariant under spatial rotations and translations x — > i?x + d, where R is an S'0(2)-matrix 
and d a translation vector in R 2 , and under global SO(2) iso rotations defined earlier ([]]). 
They are also invariant under simultaneous reflections in the Euclidean plane and on the 5| 
manifold: 

(x 1 ,x 2 ) -> (-2:1,2:2) 

(01,02,03) -> (-01,02,03) • (6) 

Physically one may think of this transformation as simultaneous electric charge conjugation 
and parity operations. 

Translationally invariant fields necessarily have degree zero, but it is possible to write 
down fields of arbitrary degree which are invariant under the reflection (|j) and under a 
combination of a rotation by some angle \ G [0, 2ir) and an iso-rotation by — n\- The 
appropriate ansatz for the scalar field 4> is, in terms of polar coordinates (r, 6) for x, 

/ sin/(r) cos n6 \ 
(f)(r,9) = sin/(r) smn9 . (7) 
V cos/(r) / 
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This is a two-dimensional version of the hedgehog ansatz used in the three dimensional 
Skyrme model. Here we consider configurations with their symmetry center, defined by 
4> = — n, at the origin. 

Under the reflection (§) the gauge field transforms as 

A) i ► -Aq, A r ^ -A r , A e <-> A g , (8) 

with the polar and radial coordinates of Ai defined analogously to those of ji before Eq. ([□]). 
Thus the requirement of rotational symmetry and reflection symmetry implies the following 
form for the gauge potential 

A = A r = , A e = na(r) , (9) 

where a is an arbitrary function and the factor n is introduced for convenience. For such a 
gauge field the electric field vanishes, and the magnetic field is given by 

B = n-. (10) 
r 

To ensure that the field is regular at the origin we impose 

/(0) = Jbr, A;eZ and o(0) = , (11) 

and the finite energy requirement implies for the function / 

hm/(r) = 0. (12) 

With these boundary conditions the topological charge Q of the hedgehog configuration (0) 
is equal to — n if k is odd and zero otherwise |J. In the following we will restrict attention 
to k = 1. 

For configurations of the form (0) and (^) the current j a has only one non-vanishing 
component, namely j$: 

jo = n(l + a)(l + A 2 /' 2 ) sin 2 /. (13) 

The electric field vanishes and the magnetic field is independent of 9, so only the ^-component 
of the Euler-Lagrange Eq. (H) is non-trivial. Thus the field equations (H) and @ imply two 
equations for a and /, which read as follows: 

f' 

f" (l + A V sin 2 /) + — ((2(5r)' - o) A 2 a sin 2 / + A 2 r/ a 2 sin / cos / + l) 

— a 2 sin / cos / — ji 2 sin / = , (14) 

where a = n(a + l)/r, and 

a" - -a = g 2 (l + a)(l + A 2 /' 2 ) sin 2 / . (15) 
r 
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Note that a = is not a solution of the second equation. Since other constant solutions are 
forbidden by the boundary condition (|II|) it follows that all solutions will have a non-trivial 
magnetic field given by fliCf). 

We will discuss the solutions of ( |14D and (|l^) in detail in the next sections, but first we 
want to address another question of principal interest. Are there finite-energy solutions of 
the field equations (§) and (H) which have a time-independent purely radial electric field? 
First we note that, in two spatial dimensions, finite energy solutions necessarily have zero 
electric charge 

q = Jd 2 xj (x). (16) 

For it follows from Gauss's law (|i~0|) that the modulus of the electric field falls off like q/r 
for large r. Hence the electric energy J d 2 x Ef diverges if q 7^ 0. However, this argument 
does not rule out finite energy solutions with a non-trivial but spherically symmetric charge 
distribution which integrates to zero. Like for example the hydrogen atom such a distribu- 
tions would only produce a short range electric field. We claim that this possibility is not 
realized in our model. To prove this assertion we must allow for more general fields than 
considered so far. In particular we can no longer impose the reflection symmetry since 
it eliminates radial electric fields from the start. Imposing invariance under simultaneous 
spatial rotations and iso-rotations leads to the following general form for the scalar field: 

' sin/(r) cos(n9 — x( r ,t)) \ 
4>(r,6,t) = sin/(r) sin(n0 - x(r, t)) , (17) 
V cos f(r) ) 

where x( r > t) is an arbitrary function of r and t. However, having introduced this function 
we can immediately remove it by a gauge transformation which brings ( [1?] ) into the standard 
hedgehog form (|7j) . Thus having fixed the gauge we write down the most general gauge field 
which gives rise to a purely radial time-independent electric field 

A = v(r), A r = h(r)t } A e = na(r), (18) 

where v and h are arbitrary functions of r. The electric field is then 

£ 4 = -K(r) +%))-• (19) 

r 

Inserting this ansatz into the field equation (H) and © leads to a complicated set of coupled 
differential equations. Let us first consider the "electromagnetic" equations @. The 6- 
component of the current j a is still given by (O), so the equation implied by the ^-component 



of (|9]) is (15) as before. However, both the t and r component of @ now lead to non-trivial 
equations which read 

thsin 2 f = (20) 

v" + -v' = g 2 v(~L + / 2 )sm 2 / • (21) 
r 
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The first clearly implies that h is identically zero. To analyze the second we first note that 
v has to satisfy the boundary condition v'(0) = to ensure that the electric field is regular 
at the origin and that for large r, v'(r) has to tend to faster than 1/r for the electric 
field energy to be finite. However, under these conditions we can multiply by (rv (r)), 
integrate both sides over r from to oo and finally integrate by parts to obtain 

J™ rdr ([v') 2 + v 2 g 2 (l + A 2 /' 2 ) sin 2 /)) = . (22) 

Since the integrand is positive definite it follows that v = everywhere. Thus v = h = 0, 
and the functions / and a satisfy the same equations as before. In particular, the electric 
field of the solution vanishes everywhere. 



IV Asymptotic Properties 

To learn more about the minimal energy solutions of the rotation and reflection symmetric 
form (0) and @ we need to solve of the boundary value problem posed by the coupled second 
order equations (|i~3[ ) and (|I5|) and the boundary conditions fllT]) and ([12]) . This requires a 
careful analysis of the equations near the regular singular points r = and r = oo of (|1~4"D 
and (0). 

At the origin, / and a behave as follows 

fir) « vr + c r n + c 2 r n+2 (1) 
a(r) « ^ + 9 ^^ r 2 ^ 2 , (2) 

where 5\ is the Kronecker symbol, Co and c? are arbitrary parameters and c 2 is a function 
of c , G? , n and /x. 

We already know that / tends to zero for large r ([12]). Thus the Eq. (|1^) becomes, for 
large r, 

ra = a . (3) 

which is solved by a constant function or by a(r) = r 2 . Since the latter leads to a magnetic 
field with infinite energy we conclude that there exist a number such that 

lim a(r) = a^. (4) 



Note that the finite energy requirement does not impose any restrictions on the value of a^. 
Since is related to the magnetic flux 

$ = [ d 2 xB (5) 
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via 

$ = 27ma OQ , (6) 

there are also no a priori restrictions on the value of the magnetic flux. This should be 
contrasted with the situation in the ablian Higgs model, for example, where the flux is 
quantized. 

It follows from (|) that the Eq. @ can be simplified for large r, and becomes 



r+j- r ( V 1) ' +p')/ = o. (7) 



The solutions of this equation are the modified Bessel functions K m (fir),m = n^a^ + 1). 
Thus / is asymptotically proportional to 



We then deduce from ([15]) the asymptotic proportionality 

fl(r)<x-i«-*-(l-o(£)). (9) 

This shows in particular that the magnetic field has no long range component. 

Having understood the asymptotic properties of Eqs. (|14D and ( |15"D it is relatively 
straightforward to solve them numerically. We have done this for n = 1 and n = 2 us- 
ing both a shooting method and a relaxation technique, with identical results. The more 
general study of static multisoliton solutions in the (ungauged) baby Skyrme model in || 
suggests that minimal energy configurations have the rotationally symmetric form consid- 
ered here for degrees 1 and 2, but are less symmetric for higher degrees. We expect that 
the rotationally symmetric solutions of degree n > 2 are similarly not true minima of the 
potential energy in our gauged model, and we therefore do not consider them here. 

The solution for n = 1 corresponds to the basic soliton of our model, which we call a 
gauged baby Skyrmion. Already the simple asymptotic analysis carried out in this section 
tells us some of its qualitative features. The profile function / of the gauged baby Skyrmion 
has the same asymptotic behaviour as that of the the ungauged baby Skyrmion discussed in 
||. As explained there the resulting asymptotic forms of the fields <f>i and 02 are the same 
as those produced by two orthogonal scalar dipoles in classical Klein-Gordon theory. Thus, 
from the point of view of the scalar (or matter) fields, the gauged baby Skyrmion looks from 
afar like a pair of orthogonal scalar dipoles in the plane of motion. 

In addition, however, the gauged baby Skyrmion has a non-trivial electric current distri- 
bution with a magnetic dipole moment orthogonal to the plane of motion. Such a magnetic 
dipole moment does not produce a long range magnetic field in (2+l)-dimensional electro- 
magnetism, but it does carry magnetic flux. Our soliton similarly has no long-range magnetic 
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field and also carries magnetic flux. Thus from afar a gauged baby Skyrmion looks like a 
triplet of mutually orthogonal dipoles: two scalar dipoles in the plane of motion and one 
magnetic dipole orthogonal to it. 



V Numerical Results 

To compute explicit solutions of fll^D and fllS]) we need to fix the parameters of the model. 
We can fix the energy and length scales by setting H = 1 and A = 1, so that we are working in 
geometric units where all quantities are dimensionless. We further choose /i 2 = 0.1 in order 
to be able to compare our results with the discussion of the ungauged baby Skyrme model 
in H (where the reader can read about the motivation for this choice). The electromagnetic 
coupling constant g, however, is not kept fixed. In fact we are particularly interested in the 
dependence of the baby Skyrmion's properties on g. Consider first two global properties: 
the energy or mass, and the magnetic flux $ fl5|). 

The energy of the solutions for n = 1 and n = 2 is shown in Fig. 1 as a function of g. 
When n = 2, the energy is less than twice the energy of the n = 1 soliton for all g, so this 
solution may be thought of as a bound state of two gauged baby Skyrmions. Recall that the 
Bogomol'nyi bound on the energy depends on whether g < \x or g > \x. In the first regime 
the bound is independent of g, but in the second it decreases like 1/g 2 , compare (f|) and 
(|5p. From Fig. 1 it is clear that the energy shows a very similar dependence on g. For both 
n = 1 and n = 2 it is essentially constant in the first regime, staying approximately 50 % 
above the Bogomol'nyi bound. Here we find in particular that in the limit g — *> the energy 
of the n — 1 and n = 2 solution tends to E\ = 1.564 ■ Att and E 2 = 2.936 • 4tt respectively, 
which agrees with the calculation for the ungauged model in [[|. In the second regime, by 
contrast, the energy, like the Bogomol'nyi bound, decreases rapidly as g is increased further. 
However, our numerical results suggest that for both n — 1 and n = 2 the energy tends to a 
non-zero limit for g — > oo. 

The precise dependence of the magnetic flux $ on the coupling constant g is shown in 
Fig. 2 for both n = 1 and n = 2. In the limit g — > the magnetic flux tends to zero, which is 
what one expects physically and which one can understand analytically by noting that in the 
limit g — > Eq. (p~5f ) becomes Eq. (|3]). As we saw in our earlier discussion of that equation, 
the only finite energy solution is the constant solution. It then follows from the boundary 
condition ( O) that a and hence also $ vanish in this limit. Furthermore, integrating Eq. 



(|i~5|) once and using lim^oo B(r) = we find 



5(0) = g 2 f dr (i±^j (1 + A 2 /' 2 ) sin 2 /. (1) 

Numerically we observe that the dependence of a and / on g is small for weak coupling, and 
that the flux $ is approximately proportional to B(0), so we expect $ to grow quadratically 
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with g for small g. This is precisely what the double logarithmic plot in Fig. 2 shows: there 
is a weak coupling range (0 < g < fi) where log$ is a linear function of logg with gradient 
2. Thus in this region the magnetic flux of the solutions of degree n = 1 and n = 2 is 
approximately 

®n-C n g\ (2) 

where C x ~ 24.5 and C 2 « 31.5. 

For large g the flux of the degree n solution tends to —27m for n = 1,2. These are 
precisely the allowed values of the magnetic flux in models such as the abelian Higgs model 
where the flux is quantized for topological reasons. Thus, although there is no such reason 
for flux quantization in our model we observe an effective quantization in the strong coupling 
limit. 

To understand the effective flux quantization and the limit of the energy as g — > oo 
better, we look at the dependence of the gauged baby Skyrmion's shape on the coupling g. 
The function a characterizing the magnetic field is plotted in Fig. 3(a) for several values 
of g. Note in particular that for strong coupling, a tends to a step function, taking the 
value at the origin but —1 everywhere else. Thus at strong coupling the magnetic field is 
increasingly localized at the origin. This is certainly consistent with Eq. (O) in the limit of 
large g, although it is not obviously implied by it. Note also that = —1 implies, via @, 
our earlier observation that at strong coupling the magnetic flux is quantized in units of 2ir. 

In Fig. 3(b) we plot the profile function / for a range of couplings. We have not included 
more plots at weak coupling (g < fi) because the profile function barely changes in this 
regime. Thus, as assumed in the calculation of the proton-neutron mass difference in Q , the 
back reaction of the electromagnetic field on the scalar field is negligible at weak coupling. 
At strong coupling, however, the profile function changes significantly, and the configuration 
becomes more localized. We conjecture that / also tends to a singular step function in the 
limit g — > oo. To justify this conjecture, consider Eqs. flUD and dl5|). When a is the step 
function described above, these equations decouple everywhere except at the origin, and the 
first equation becomes the Euler-Lagrange equation derived from the functional 

F[f] = Jdrr(±f a + ii 2 (l-coBf)y (3) 

However, this functional cannot be minimized by a non-singular function satisfying the 
boundary conditions (|H|) and (12). For suppose there were such a minimum, and call it f\. 
Then define f K {r) = and find that F[f K ] < F[fi] for k < 1. Thus one can always 

lower the value of F[f] by making / spikier. 

To sum up, we have the following picture for the n = 1 and the n = 2 solitons in the 
strong coupling limit: both the magnetic field and the energy distribution become localized 
near the origin, tending to singular distributions as g — > oo. In this limit, the total energy 
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does not vanish because of the contributions from gradient terms in the energy density, and 
the magnetic flux tends to —27m. 

To end, let us look at the solutions of ([T3|) and (|T5) for a particular value of g in more 



detail. Ideally we would like to pick a "physical" value for g which allows us to estimate 
the magnitude of the effects we have observed in the physically relevant (3+l)-dimensional 
Skyrme model. In practice it is not clear how to define "physical" here, but the following 
procedure should lead to a value for g which is at least in the right ball park. 

The basic idea is to treat our model as if its solitons described physical baryons and its 
elementary quanta were physical pions: this allows us to fix the energy and length scale, and 
to compute a definite value for the coupling constant g from the physical value of the fine 
structure constant. Thus we assume, in accordance with nature, that the electromagnetic 
coupling constant is small. 

First we fix the energy scale H by identifying the mass of a gauged baby Skyrmion 
with the physical nucleon mass of 940 MeV. Since for small g the mass of the gauged baby 
Skyrmion is virtually independent of g, we pick the value of the mass at g = 0; this leads 
to H = 48 MeV. To find a physical length scale A we note that 1/fi is the equivalent to 
the Compton wavelength of the pion in the Skyrme model. Thus we choose A such that 
= 1.41 MeV, i.e. A = y/0.1 ■ 1.41fm = 0.45 fm. To compute Planck's constant in 
geometric units we write 

h = 197.3 MeV fm « 9.1(48MeV)(0.45 fm) (4) 

and deduce that H = 9.1 in geometric units. Finally we use the physical value of the fine 
structure constant a = e 2 /(4nh) »s 1/137. Here e is the electron's charge which is related to 
the coupling constant g via e = gh. With h = 9.1 we conclude g = 0.1 in geometric units. 

At this value for g the gauged baby Skyrmion is lighter than the baby Skyrmion at 
g = by AE\ = 0.12, which is 5.9 MeV in physical units. For the solution with n = 2 the 
corresponding energy difference is AE 2 = 0.56, which is 27 MeV in physical units. It is also 
interesting to look at the g-dependence of the difference between the energy of the n = 2 
solution and twice the energy of the n = 1 solution, which may be interpreted as a binding 
energy. In the g = case, this is about 6.6%, but it is 14% when g = 0.1. Thus the inclusion 
of the electromagnetic field leads to a more strongly bound n = 2 soliton. 

The precise shape of the solutions for g = 0.1 can be seen in Fig. 4., where we plot the 
energy density and the magnetic field for both n = 1 and n = 2. The baby Skyrmion's 
energy distribution is bell-shaped and peaked at the origin, whereas the energy distribution 
for the soliton of degree 2 is maximal on a ring with radius r = 1.78. 
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VI Conclusions 



In this paper we have studied soliton solutions of the coupled Skyrme-Maxwell system in 
(2+1) dimensions. The rotationally symmetric solitons we have considered necessarily carry 
a magnetic field but the electric field is zero. The magnetic flux can take arbitrary values, 
but in the strong coupling limit we observe an effective flux quantization. The soliton 
mass decreases when the electromagnetic coupling constant is increased and all the other 
parameters of the model are kept fixed. Thus a baby Skyrmion can lower its mass by 
interacting with the electromagnetic field. 

Although the U(l) gauge group is unbroken the baby Skyrmions' magnetic field is short- 
ranged. The reason for this is that the electromagnetic current carried by the baby Skyrmion 
only has a magnetic dipole component; in (2+1) dimensions static magnetic dipoles, however, 
have no long-range fields in Maxwellian electromagnetism. This observation has important 
consequences for the interaction of gauged baby Skyrmions. Since the scalar fields fall off 
like exp(— fir) and the magnetic field like exp(— 2/xr), the magnetic forces will be negligible 
compared to the scalar forces between two well-separated gauged baby Skyrmions. Thus, 
to first approximation, the forces should be the scalar dipole-dipole forces between purely 
scalar baby Skyrmions discussed in detail in ||. 
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Figure Captions 



Fig. 1: Energy "per Skyrmion" as a function of the coupling constant g for n — 1 and n = 2 
. The dashed line is a plot of the Bogomol'nyi bound, see (f|) and (^). 

Fig. 2: The magnetic flux "per Skyrmion" as a function of the coupling constant g for 
n = 1 and n = 2. 

Fig. 3: The function a (a) and the profile function / (b) of the n — 1 solution for various 
values of the coupling constant g. 

Fig. 4: The energy density (a) and the magnetic field (b) for n = 1 and n — 2 at g — 0.1; 
the function e plotted in (a) is the integrand of ( |TB| ) divided by 47r. 
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